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Abstract 

It is possible to derive a monopole action from vacuum configurations ob- 
tained in Monte-Carlo simulations extending the method developed by Swend- 
sen. We apply the method to compact QED both in the Villain and in the 
Wilson forms. The action of the natural monopoles in the Villain case is 
in fairly good agreement with that derived by the exact dual transforma- 
tion. Comparing the monopole actions, we find (1) the DeGrand-Toussaint 
monopole definition may be useful for f3v larger than about 0.5, (2) the Vil- 
lain model well approximates the Wilson one for /? smaller than jSc and (3) 
in the Wilson action the monopole condensation occurs in the confinement 
phase and f3c may be explained by the energy-entropy balance of monopole 
loops like in the Villain case. 
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To understand confinement mechanism is very important but still unresolved problem 
in quantum chromodynamics (QCD). A promising idea is that the (dual) Meissner effect 
due to condensation of some magnetic quantity is the color confinement mechanism in QCD 
This picture is realized in the confinement phase of lattice compact QED 

The QED partition function in the Villain form []Tn|,|TT| is 

Zv = \{ de{s,f,) n E ^M-^ E [(^A^) - 2vrn,.(.)]^), (1) 

where the plaquette variable 6^j,{s) G (— 47r, 47r) is given by the link angles 6{s, //) G [— vr, vr) 
as 9^y{s) = 9{s, fi) + 9{s + fi,!/) — 6{s + z>, /i) — 6{s, v). A dual transformation can be done 
exactly , leading us to a partition function describing a monopole Coulomb gas | 



Zm = {X{ E )(n'^5;,™,(s),o)exp(-27r2/3v^ ^ m^{s)D{s - s')m,{s')), (2) 

m^(s)=— oo ^ s,s',ij, 

where D{s — s') is the lattice Coulomb propagator and d' is the backward difference. Here 
mfj,{s) is a monopole current: 

m^is) = ^e^^p„d^np^{s + /i), (3) 

which is defined on a link of the dual lattice (the lattice with origin shifted by half a lattice 
distance all in four directions), d is the forward difference. Note that the coupling constant 
in (^ satisfies the Dirac quantization condition between electric and magnetic charges. One 
can prove the monopole condensation from energy-entropy balance using the action (|^) and 
the (bulk) entropy ~ Lxln7 where L is the monopole loop length [^^-j^. The critical /3y 
is estimated to be py = ln7/(27r^D(0)) = .637 0] which is in good agreement with the 
Monte-Carlo data /9{> = .645 H . 



On the other hand, the Wilson partition function of lattice compact QED is 

Zw = l[ r d9{s, fi) expi-Pw E (1 " cos9^,is))]. (4) 

Since the action of the Wilson form is the simplest and is used commonly also in non-abelian 
gauge theories, it is important to study the confinement mechanism of compact QED with 



this action. However the exact dual transformation which enables us to make the above 
analytic study is impossible in this case, although this action (^) is easy to perform Monte- 
Carlo simulations. To speak rigorously, hence, the monopole condensation is not yet proved 
in the case of the Wilson action. 

Using the Fourier transformation of (^) and the property of the modified Bessel function, 
the Wilson action is approximated for very large (3 and for small f3 by the Villain form (P. 
How good the Villain approximation to the Wilson form is was studied in |llO| , p!I| extensively. 



In Ref. [0, physical quantities like internal energy and specific heats can be well reproduced 
for small f3v up to f3y, if both the coupling constants satisfy 

where Jq and Ii are usual modified Bessel functions. The approximation becomes rapidly 
worse above the transition point except for very large (3. For moderately large f3v, the Villain 
partition function gives results rather similar to a mixed cosine action with a modified 
coupling relation. However even if some physical quantities are well reproduced by the 
Villain approximation, it does not always mean that the confinement mechanisms of both 
models are the same. Really it is pointed in that there are physical quantities not well 



reproduced even for small (3. Considering applications to QCD, we want to know directly 
the confinement mechanism in the Wilson action. 

For that purpose, one has to find first monopole currents in the Wilson action where there 
is no natural definition like contrary to the Villain case. The only known definition is 
that given by DeGrand and Toussaint 0. The plaquette variable 9^^{s) can be decomposed 
into 

Of^u{s) = 6f,u{s) + 2Tcm^^{s), (6) 

where 9^u{s) G [— 7r,7r) and m^i,{s) can be regarded as a number of the Dirac string pene- 
trating the plaquette. DeGrand and Toussaint pi defined a monopole current as 
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^^Ji [s) = ^^fj,iypadump„{s + (j). (7) 



Recently Schram and Teper studied how well the DGT monopoles m^*-^(s) approxi- 
mate the natural ones m^{s) based on Monte-Carlo simulations of the Villain form. They 
suggested the DGT prescription can be meaningfully used in the studies of the U(l) phase 
transition. 

Next we have to find a partition function described in terms of the monopole currents 
corresponding to the dual form of the Wilson action in order to study the energy entropy 
balance. Actually an exact dual transformation is impossible in almost all models where 
dual variables are expected to play an important role in the dynamics. It is very important 
to develop a method determining a dual theory for such models. It is the most important 
original point of this note to give a method fixing a monopole action from a given ensemble 
of configurations generated in Monte-Carlo simulations. To demonstrate the usefulness of 
our method, we first apply the method to the natural monopoles in the Villain case of U{1) 
lattice gauge theory where the exact monopole partition function is known as in (0). Next 
we derive the action for the DGT monopoles in the Villain case and test reliability of the 
DGT prescription. Finally we apply the method to the Wilson case to check the occurence of 
the monopole condensation. A brief description of our method and preliminary applications 
were given also in |]n|,|T3|. 



We extend the method of Swendsen fl^ to determine an action from a given ensemble 
of monopole loops in vacuum configurations. A theory of monopole loops is given in general 
by the following partition function 

oo 

Z = C[l E )iI[^0',k,is),o)eM-S[k]), (8) 

^'f" fc;i(s)=— oo ^ 

where fc^(s) is the conserved integer-valued monopole current defined above in (^) or (|^ and 
S[k] is a monopole action describing the theory. Consider a set of all independent operators 
which are summed up over the whole lattice. We denote each operator as Si[k]. Then the 
action can be written as a linear combination of these operators: 

s[k] = j2mk], (9) 



where fi are coupling constants. Really the monopole partition function takes this form. 
The expectation value of an operator of monopole currents are estimated as 



Let us now determine the monopole action using the monopole current ensemble which 
are calculated by vacuum configurations generated by Monte-Carlo simulations. Since the 
dynamical variables here are kf^{s) satisfying the conservation rule, it is necessary to extend 
the original Swendsen method. Consider a plaquette (s', /i', u') instead of a link on the dual 
lattice. Define Si[k] as a part of an operator Si[k] containing the monopole currents on the 
links in the plaquette chosen, i.e., k^i{s'), k^i^s' + /i'), k^i{s' + v') and kyi^s'). Then we get 

oo 

(n E ){X{6d',k,is)fl)Sa[k]eM-iimk]) 

s,M fc^(s)=-oo s i 

I oo I oo oo oo oo 

= (n E )(n^a;..w,o)[ E E E E 

s,iJ. k^{s)=-oo s k^,{s')=-oo k^,(s'+il')=-ook,(s'+u')=-ook^,{s')=-oo 



i 

exp{-J2MSi[k]-S,[k])), (11) 

i 

where the primed product means removing the sites and the links in the plaquette chosen. 
Note here that the sum of the current conservations on the four sites d'^k^{s') + d'^k^{s' + 
/^') + '^'fj.kfj.is' + ly') + d'^k^{s' + /i' + v') does not contain any current on the four links of the 
plaquette adopted. Hence 

/ oo / 

(11) = (IT E )([\^d[^kt^{s)fi)^dlk^,(s')+d'^k^,{s'+ii')+d'^k^,{s'+C')+d'^k^,{^ 
•'iM fcp(s)=— oo ^ 

exp(- ^ U{Sm - ^.[^]))(E ^)kSa[k] exp(- ^ fA[k]) (12) 

i i 
oo 

= (n E msd'k,is),o)m^M-j:mk]), (13) 



^lA* kfj,{s)=—oo ^ 



where 



g r^. ^ {J:5)i,Sa[k,{ky]expi-E^mk{ky]) 
{ES)kexp{-J:^mk,{ky]) 



Here {k}' do not contain the four currents on the hnks of the plaquette considered and 

oo oo oo oo 

(E%-_ E ^ E _ E ^ E 

k^i{s')=—oo fc^; (s'+/i')=— CO fc^/ («'+!/')=— CO fe^/(s')=— oo 



Since there are current conservations at all sites in (|l^), we get 

d'Ms) = d'^k,{s') + k,is') + his') - k,is') - k,,{s') (16) 

= k^i{s) + kui{s') - k^i{s) - ky>{s). (17) 

Similarly, we see 

&Jz^{s' + ii') = his' + ii') - hi^') - ku'is' + fi') + V(s'), (18) 

d'^hi^') + d'Ji^is' + ?) = h is' + z^O + V (^0 - h is' + - h is) ■ (19) 

Now use the following relation: 

CO 

E h^,{s%k^,{s')+M = 1- (20) 

M=-oo 

Then 

oo oo 

iT.^)kSa[kis'),K'is' + ^i'),h,is' + l)'),his'),{k}'] = Y: E ^V(^').V(^')+M 

M=-oofc^,(s/)=_oo 

oo oo oo 

^ E ^k^>{s'+il'),k^,{s'+il')+M E ^k^^,is'+u'),k^,{s'+C')-M E h^,{s'),k^,{s')-M 

k^i{s'+^')=—oo k^i {s'+u')=—oo fc^/(s')=— oo 

X Sa[k^'is') + M, k,,is' + [i') + M, V(s' + ?) - M, K,is') - M, {A;}'] (21) 

CO 

= E Sa[k^,is') + M,k,\s' + ii') + M,k^,is' + u')-M,k,,is')-M,{k}']. (22) 

Af=-oo 

When the DGT monopoles are used, the sum with respect to M is restricted from the 
minimum mi to the maximum m2, where 

mi = -2 - Min{ V (s') , (s' + h) , " V (§' + u'),-Kis')}, (23) 
m2 = 2 - Max{ V (s') , k,. (s' + fi'), -V (s' + /^') , -K' (s') } • (24) 



Note that the DGT monopoles take integer values between -2 and 2 only as seen from the 
definition. Then using (p!0|), we get 

{Sa[k]) = {Sa[k]), (25) 

where 

g iji ^ Sm=-oo '5'a[fe] exp(— fiSi[k]) 

and 

k^{s) = k^{s) + M{6s,s'Sf,,f^' + S^ ,.,^^,6^y - 5^ .,,+^-,5^,^' - Ss,s'Si,y). (27) 

Introducing a new set of coupling constants {/«}, we define 

^ ^ EM=.ooS^[k]exp{-E.mk]) .28) 
EM=-ooexp(-Ej.^.[^]) ' 



where k is defined in (^). When all /j are equal to /j, one can prove an equality from (^51) 

{S^k]) = {S^k]). 

When there are some fi not equal to fi, one may expand the difference as follows: 

{Sa - Sa) = Y^iS^S, - SaS,){h - /fc), (29) 
b 

where only the first order terms up to 0{fb — fb) are written down. This allows an iteration 
scheme for determination of the unknown constants fi from the ensemble of {A;^(s)}, which 
are generated in Monte-Carlo simulations. 

Practically we have to restrict the number of interaction terms Q. We adopted quadratic 
interactions of up to 32 types listed in Table | in these studies. The first six terms are shown 
also graphically in Fig. |l|. 



^ All possible types of interactions are not independent, since d'^k^{s) = 0. We can get rid of 
almost all interactions between different components of the currents from the quadratic action by 
use of the conservation rule. 



Before getting the configurations, we checked the auto-correlation time and thermahza- 
tion. For thermahzation, we need more than 500 iterations in the Villain case, whereas 50 
iterations are enough in the Wilson case near the critical (3. The auto-correlation quickly 
disappears in the Wilson case, whereas the situations in the Villain case are much worse 
near the critical (3. Hence using the Villain action, we generated 100 gauge field configu- 
rations separated by 50 sweeps after a thermahzation of 1000 sweeps for (3v < 0.56. At 

= 0.62, we separated 300 iterations. For more than the critical (3y ~ 0.64, we performed 
100 independent runs with different initial values and adopted one configuration after 1000 
thermahzation loops for each run. Each run in the Villain case was done on 8^ lattice. In the 
Wilson case, we took 100 configurations separated by 50 sweeps after 1000 thermahzation 
sweeps for every on 8^ and 12^ lattices. The monopole currents are defined following @ 
and (0). The statistical errors were estimated with the jackknife method. 

Our results are the following. 

1. The first coupling constants /i ~ /g of the action for the natural monopoles in the 
Villain case are plotted in Fig. ^ in comprison with the theoretical values given by (^. 
/i agree well with those of the theoretical values, whereas there are small discrepancies 
with respect to /2 ~ /e- The discrepancies may come from the truncation of the terms 
of the action taken. At the critical /Jy ~ 0.64, fi crosses the ln7 entropy line. 

2. The same coupling constants of the action for the DGT monopoles in the Villain case 
are plotted in Fig. ^ Both /i of the natural and the DGT monopole actions are 
compared in Fig. |. We see there is a large deviation for small (3. For jSy > 0.5, we 
may use the DGT prescription. Schram and Teper [|1^] recently showed the difference 



of the natural and the DGT monopoles can be reproduced by a random distribution 
of dipoles having trivial long-distance behaviors. They concluded therefore the DGT 
prescription can be used for /Sy > 0.3. /i is the dominant part of the action and 
it plays an important role in the energy-entropy balance. There is a large gap with 
respect to /i at f3y = 0.3. It seems dangerous to use the DGT definition for such small 
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I3v- Fortunately the DGT prescription looks rather good around the ctitical I3y. 

3. Using the DGT definition, we get the monopole action in the Wilson case. The first 
six coupling constants of the action on 12^ lattice are plotted in Fig. ^. There is a 
small volume dependence for (3\y larger than the critical value( ~ 1.0). /i on 8^ and 
12^ lattices are compared in Fig. |^. The energy of a monopole loop of length L may 
be approximated by a self-energy part fiL when L is large P,p^. The entropy line 
ln7 and the /i line crosses at (3w = 1-05 which is very near to the critical Pw ~ 1-0. 
Namely the critical coupling constant determined from the monopole condensation due 
to the energy-entropy balance of monopole loops agrees to that of the deconfinement 
transition. The monopole condensation really occurs in the confinement phase also in 
the Wilson model. This is the first direct evidence for the occurence of the monopole 
condensation in the Wilson action without the use of the Villain approximation. 

4. Similarity of the Villain and the Wilson partition functions is studied by comparing 
both /i in Fig. 0, where Pw is transformed into Pv using (||). Both are in good 
agreement for small (3 less than the critical value ~ 0.64. For larger /3, both begin to 



deviate rapidly. This is consistent with the conclusion given in |11 



Finally we have seen that our method is very useful in the study of compact QED. The 
method can be applied also to more interesting QCD cases. For preliminary reports, see 



l^ , p!5[| . The detail of the results will be published elsewhere |]T8 . 

We wish to acknowledge Yoshimi Matsubara for useful discussions. This work is finan- 
cially supported by JSPS Grant-in Aid for Scientific Research (B) (No. 06452028). 
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TABLES 

TABLE I. The quadratic terms of the monopole action adopted. Only the partner of the 
current multiplied by k^{s) are listed. All terms in which the relation of the two currents is 
equivalent are added to make each Si [k] invariant under translation and rotation. Here d = z> + p, 
b = !> + p + u, c = fi + C' + p + u, d = 2/1 + i>, e = 2z) + /x and / = 2z) + a), where jl, 0, p and a) 
denote unit vectors in four different directions. 
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FIGURES 



FIG. 1. The first six terms of the monopole action adopted. 

FIG. 2. Couphng constants /j versus (3 in the Villain model when the natural monopole 
currents 'm^{s) are used. The solid (dotted, dashed and dott-dashed ) line denotes the theoretical 
curve of /i(/2 = /a, Ja = h and /e). 

FIG. 3. Coupling constants fi versus [3 in the Villain model when the Degrand-Toussaint 
monopole currents m^'^{s) are used. The lines show the theoretical curves predicted by the 
monopole action for the natural monopoles. 

FIG. 4. Coupling constants f\ versus (3 in the Villain model for both the natural and the 
Degrand-Toussaint monopole currents. 

FIG. 5. Coupling constants fi versus (3 in the Wilson model on 12^ lattice where the De- 
grand-Toussaint monopole currents m^^{s) are used. 

FIG. 6. Coupling constants /i versus (3 in the Wilson model both on 8^ and 12^ lattices. 

FIG. 7. Coupling constants /i versus (3 in the Villain and the Wilson models, where the DGT 
monopoles are used. The solid line denotes the theoretical value of /i of the action for the natural 
monopoles. 
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